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ABSTRACT 

We derive the total cold gas, atomic hydrogen, and molecular gas masses of approxi¬ 
mately 24 000 galaxies covering four decades in stellar mass at redshifts 0.5 < z < 3.0, 
taken from the CANDELS survey. Our inferences are based on the inversion of a 
molecular hydrogen based star formation law, coupled with a prescription to separate 
atomic and molecular gas. We find that: 1) there is an increasing trend between the 
inferred cold gas (Hi and H 2 ), Hi, and H 2 mass and the stellar mass of galaxies down 
to stellar masses of 10^ Mq already in place at z = 3; 2) the molecular fractions of 
cold gas increase with increasing stellar mass and look-back time; 3) there is hardly 
any evolution in the mean Hi content of galaxies at fixed stellar mass; 4) the cold gas 
fraction and relative amount of molecular hydrogen in galaxies decrease at a relatively 
constant rate with time, independent of stellar mass; 5) there is a large population of 
low-stellar mass galaxies dominated by atomic gas. These galaxies are very gas rich, 
but only a minor fraction of their gas is molecular; 6) the ratio between star-formation 
rate (SFR) and inferred total cold gas mass (Hi -h H 2 ) of galaxies (i.e., star-formation 
efficiency; SEE) increases with star-formation at fixed stellar masses. Due to its sim¬ 
plicity, the presented approach is valuable to assess the impact of selection biases on 
small samples of directly-observed gas masses and to extend scaling relations down 
to stellar mass ranges and redshifts that are currently difficult to probe with direct 
measurements of gas content. 
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1 INTRODUCTION 

Observations in the local Universe have revealed that star 
formation (SF) is closely linked to the density of molecular 
gas. SF in the Milky Way takes place in dense, massive and 


cold giant molecular clouds (Solomon et al. 1987 McKee 


& Ostriker 2007[ Bolatto et al. 2008). Recent works have 


emphasized that there is a strong correlation between the 
star-formation rate (SFR) density and the density of molec¬ 
ular hydrogen (H 2 ), while the correlation with the density 
of atomic hydrogen (Hi) is weak or absent (Wong & Blitz 


2002 Bigiel et al.||2008[ |2011| Schruba et al. 2011). These 

results have stimulated the idea that not all the cold gas 
in a galaxy is necessarily available for SF. A proper under¬ 
standing of the evolution of the atomic and molecular gas 
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content of galaxies is a key ingredient that will shed light on 
the physics that regulates the stellar buildup of galaxies. 

Current surveys of the cold gas content of galaxies at 
high redshift only probe the molecular gas mass and are usu¬ 


ally limited to massive galaxies with high SFR (e.g., Daddi 


et al.|201^|Tacconi et al.|2010|jGeach et al.|2011||Bauernieis- 

ter et al.||2013| [Tacconi et al.pOl^ Santini et al.||2013[ ). ft 

is crucial to include the contribution from atomic hydrogen 
to thoroughly understand how the total cold gas properties 
of galaxies regulate the SFR. Furthermore, for a complete 
assessment it is necessary to study a galaxy sample that is 
not biased to the most actively star-forming objects. It is 
hoped that facilities such as ALMA (Atacama Large Mil¬ 
limeter Array), NOEMA (NOrthern Extended Millimeter 
Array), SKA (Square Kilometer Array), MeerKat (Karoo 
Array Telescope) and ASKAP (Australian SKA Pathfinder) 
will reveal the HI and H 2 content of representative samples 
of high-redshift galaxies. 

Theorists have made considerable progress develop¬ 
ing models that track the HI and H 2 content of galaxies 
pbreschkow Rawlings |2QQ9| |Eu et al.| 201Q[ Krumholz 


De^ 2012a[ |Lagos et al.||20l'l[~ Christensen et al. 2012|^ 


Kuhlen et al.||201^ [Dave et al.||2013| [Popping, Somerville| 

Trager||2014| ^inerville. Popping Trager||2015). These 


models have proven successful at reproducing the available 
observational estimates of the overall HI and H 2 properties 
of local and high-redshift galaxies. Nevertheless, observa¬ 
tional constraints at high redshift are still very limited and 
do not probe the wide parameter space covered by the mod¬ 
els. Additional information on the gas content of galaxies 
will be crucial to break the degeneracies in different phys¬ 
ical mechanisms that are included in models (e.g., SE and 
stellar feedback). 

In the meantime, we can obtain indirect constraints on 
the gas content of galaxies by using the empirical relation 
between SER density and gas density. This approach has 
been used extensively by inverting the Schmidt-Kennicutt 
relation ( Kennicutt| 1998b hereafter the KS-relation), which 


relates the SER surface density to the combined atomic and 


molecular hydrogen surface density (e.g.. 

Erb et al. 

2006 

Mannucci et al.[[2009 Troncoso et al.[[2013 

). Popping et al. 

(2012) were the first to use an inverted molecular-gas-based 


SE law in combination with a recipe to separate atomic from 
molecular hydrogen ( Blitz fc Rosolowsky||2006 ) to estimate 
the total cold gas and molecular gas content of galaxies. This 
approach was motivated by observations demonstrating that 
SER surface densities correlate almost linearly with molec¬ 
ular gas surface density (even in the low gas surface density 
regime), whereas the KS-relation breaks down at low To¬ 
tal gas’ surface densiti es ([Bigiel et al.|[2QQ8 Schruba et al. 
|2011| . [Popping et al. (2012) showed that, when inferring 
gas masses, a molecular-gas-based SE law in combination 
with a prescription to separate the atomic and molecular 
hydrogen content of galaxies is better suited to reproduce 
directly-observed gas masses and gas surface densities from 


a sample of galaxies taken from Leroy et al. (2008) than the 
total-gas KS-relation. 


Popping et al. (2012) confirmed previously observed 


trends between galaxy gas fraction and molecular gas frac¬ 
tion with stellar mass (e.g. Tacconi et al. 2010| Saintonge 


et al. 2011), through a detailed study of the inferred gas 


vey Scoville et al.||2007 ) at 0.5 < 2 : < 2.0. This initial study 
suggested that massive galaxies have lower gas fractions at 
higher redshift than less-massive objects and have lower frac¬ 
tions of their gas in molecular form. 

In this paper, we apply the met hod developed in [Pop¬ 


ping et (|2012| and updated in Popping, Behroozi & 


Peeples ( 2015| to a galaxy sample drawn from the CAN- 


DELS survey (Cosmic Assembly Near-infrared Deep Extra- 
galactic Legacy Survey; [Crogin et al.|201l|[Koekemoer et al.[ 
[2011[ ) at redshifts 0.5 < 2 ; < 3. We focus on the cold gas (HI 
+ H 2 ), Hi, and H 2 properties of galaxies over cosmic time 
and how they are related to the SE and other global prop¬ 
erties of galaxies. The CANDELS survey is deeper than the 
COSMOS survey, which allows us to study much fainter ob¬ 
jects. As such, we can probe the gas properties of the bulk 
of star-forming galaxies between z = 0.5 and 2 ; = 3. This 
cosmic epoch marks the peak in star-formation activity of 
our Universe ( Hopkins Beacom|2006 Madau & Dickinson 


2014), when the bulk of mass in today’s massive galaxies 


was formed. CANDELS provides exquisite imaging cover¬ 
ing a wide range of wavelengths to derive stellar masses 
and SER and provides reliable morphological information 
at these redshifts. The methodology enables us to infer the 
gas, HI, and H 2 masses for a large number of galaxies cover¬ 
ing a wide range in stellar masses, SERs, sizes, and redshift. 
This makes the applied method very helpful in assessing the 
impact of selection biases on much smaller samples for which 
direct gas measurements have been obtained. Eurthermore, 
it can extend scaling relations to a stellar mass and red¬ 
shift range difficult to reach through direct observations of 
gas masses. The inferred gas masses presented in this work 
have a great predictive power for future HI surveys such 
as LADUMA ( Holwerda, Blyth &: Baker[|2012 ) with instru¬ 
ments like MeerKat, ASKAP, and the SKA and future sur¬ 
veys of the molecular hydrogen content of galaxies through 
CO or sub-mm continuum imaging with for example ALMA 
and NOEMA. 

This paper is organised as follows. In Sectionwe sum¬ 
marise our method to indirectly measure the cold gas and 
H 2 content of galaxies and we present the galaxy sample 
selection from CANDELS. In Section we present our re¬ 
sults. The assumptions that were made and the applicabil¬ 
ity of our model to the CANDELS galaxy sample are dis¬ 
cussed in Section [ 4 ] We discuss our results in Section We 
summarise our findings in SectionThroughout the paper, 
we assume a A Cold Dark Matter (ACDM) cosmology with 
Hq — 70kms~^ M pc~^, ^matter = 0.28 and Qa = 0.72 ( Ko¬ 
matsu et al.[[2009' ). We assume a universal Chabrier stellar 


initial mass function (IME: Chabrier|[2003 ) and where nec¬ 
essary convert observational quantities used to a Chabrier 
IME. All reported cold, HI, and molecular gas masses in¬ 
clude a correction of 1.36 to account for helium. 


content of galaxies in COSMOS (the Cosmic Evolution Sur- 


2 MODEL DATA 

In this section we describe our method to indirectly estimate 
the cold gas (HI +H 2 ), Hi, and H 2 content of galaxies and 
the observational data to which we apply our method. 
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2.1 Obtaining indirect gas measures 

We infer the cold gas (HI + H 2 ), HI and H 2 content of 
galaxies using a combination of an empirical molecular SF 
law (based on Bigiel et al.|2008| ) and a prescription to calcu¬ 
late the H 2 fraction of cold gas ( Blitz &: Rosolowsky ||20Q6 ). 
To infer the HI and H 2 masses of a galaxy we only use the 
galaxy stellar mass, SFR, and size as input parameters. We 
pick a gas mass, distribute the gas as explained in the next 
paragraph, and then calculate a SFR for that gas mass fol¬ 
lowing the equations given below. We repeat this process 
while iterating through gas masses till convergence with the 
observed SFRs is reached. When convergence is reached we 
separate the cold gas mass into an atomic and molecular 
component using equation 

We assume the galaxy stellar mass to be distributed 
following a Sersic profile EHc(r) oc exp[—(r/r*)^/’^], where 
is the scale radius of the stellar disc for a Sersic profile and n 
the Sersic index of the galaxy. We also assume that the gas 
in galaxies is distributed following the same Sersic profile, 
with a scale radius 


T gas — Xgas T * , 


( 1 ) 


where Xgas is the scale radius of the gas disc relative to the 
stellar disc. We take Xgas = 1.7 based on a fit to the galaxy 
disc profiles presented in Leroy et al. (2008). 


We use a slightly adapted version of the star forma¬ 


tion law deduced by Bigiel et al. (2008) to allow for higher 


star formation efficiencies in high gas surface density regions. 


This is based on the results of Daddi et al. (2010) and Genzel 


et al. (2010), who found the star-format ion at high surface 


densities to follow the KS relation (a power-law slope of 1.4 
versus 1.0 for Bigiel et al. 2008). The resulting equation is 
given by 


FIsfr = 


^SF 


lOM© pc“ 


1 + 


F!h2 \ 


VSF 




/h 2 ^g 


( 2 ) 


where Esfr and Egas are the star formation and cold gas sur¬ 
face densities in M© yr“^ kpc“^ and M© pc“^, respectively; 
Aqf is the normalization of the power law in M© yr“^ kpc“^; 
Ecrit is a critical surface density; iVsF is an index which sets 
the efficiency; and /h 2 = Eh 2 /(Ehi + EH 2 ) is the molec¬ 
ular gas fraction. Following [Popping et ah (2012) we use 
Nsf = 0.5 and we take Ecrit = 70 M© pc“ . 

We use a pressure-regulated recipe to determine the 
molecular fraction of the cold gas, based on the work by 


Blitz & Rosolowsky (2006). They found a power-law rela¬ 


tion between the mid-plane pressure acting on a galaxy disc 
and the ratio between molecular and atomic hydrogen, i.e.. 


Rk: 


^ f 5^2 A ^ f ^ 
~ V^Hi J ~ \ Po 


(3) 


where Pq is the external pressure in the interstellar medium 
where the molecular fraction is unity; a is the power-law 
index; and Pm is the mid-plane pressure acting on the galaxy 
disc. We adopted Pq = 3.25 x 10“^^ergcm“^ and a = 0.8 


from Leroy et al. (2008). The mid-plane pressure can be 
described by (Elmegreen 1989| 


Pmir) — 2 GEgas(^) [Egas(^) + /(T(r)E*(r)] (4) 

where G is the gravitational constant, r is the radius from 
the galaxy centre, and fa{r) is the ratio between crgas(r) 


and crHc(r), the gas and stellar vertical velocity dispersion, 
respectively. Following Fu et al. ( |2010| , we adopt /ct(^) = 
0.1y/E*,o/E*, where E*,o = m^l{27vri), based on empirical 
scalings for nearby disc galaxies. Putting this together, we 
have an expression for the star formation surface density 
in Equation § depending on the cold gas surface density 
and stellar mass surface density. We integrate the SERs in 
the individual annuli to obtain a total SER which can be 
compared to the observed SER. 

We calibrated our method using direct observations of 
the HI and/or H 2 content of galaxies in the local and high- 
redshift Universe from [Leroy et ah] (|20Q8|), [Daddi et al. 

(|2010l. 


Tacconi et al 


Tacconi et al. 


'(2013113 


(2010), and 

Using x^-minimization, we find the best agreement between 
predicted and observed gas masses when adopting a value 
of Asf = 9.6 X 10“^ M©yr“^ kpc“^. This value is close to 


norma lizations found in Bigiel et al. (2008) and Bigiel et al. 
(2011), 8.42 X 10“^ and 4.6 x 10“^ M© yr“^ kpc“ , respec¬ 
tively. We integrate the gas disc out to 10 times rgas- We 
discuss the uncertainties on some of the individual compo¬ 
nents of this model in Section The typical systematic un¬ 
certainty of our method is 0.3 dex. 

A schematic picture of the dependencies of our model is 
presented in Eigure[^ In this figure we explore the effects of 
changing one of the three input properties of a galaxy (SER, 
size, and stellar mass) on the estimated cold gas mass and 
H 2 mass of a galaxy and the cold gas fraction (/gas) and 
molecular hydrogen fraction of the cold gas (/h 2 )- We keep 
two of the parameters fixed and vary the third. Ghanging 
the stellar mass of the galaxy has a negligible effect on the 
inferred cold gas and H 2 mass. It is only for the most mas¬ 
sive galaxies (M* 10^^ M©) that, due to the increased 

pressure from stars, the inferred cold gas mass decreases 
and /h 2 increases. Because the inferred cold gas mass re¬ 
mains relatively constant, the cold gas fraction /gas rapidly 
decreases with increasing stellar mass. 

We find a strong increase in cold gas mass, H 2 mass, 
/gas, and /h 2 when increasing the SER. This is not surpris¬ 
ing, as more gas is needed to support the higher SER. The 
increased gas mass enhances the cold gas surface density and 
pressure acting on the gas, which leads to an increase in /h 2 . 


The inferred gas masses increase as a function of scale 
radius. When increasing the scale radius, the gas densities 
decrease, therefore the molecular fraction of the gas de¬ 
creases and the SER surface density decreases as well. To 
ensure that convergence with the input SER is obtained, 
more gas has to be added to the galaxy to make up for the 
lowered H 2 fractions. This results in high gas fractions and 
low molecular hydrogen fractions. 


2.2 Data 


We infer the galaxy gas content of two deep samples of 
E160W (iL-band) selected galaxies from the GANDELS sur¬ 


vey at 0.5 < z < 3 (Grogin et al. 2011 Koekemoer et al. 


^ I Leroy et al.| (|2008HTacconi et al.| (|2010|), and |Tacconi et al.| 
( |20I3| ) assume a CO- to-H 2 conversion fa ctor of Xco = 2 x 
10^^crn“^/(K km s“^). [Daddi et al.| ( |20I0| assumes a CO-to-H 2 
conversion factor of Xqq = 2.25 X I0^^cm“^/(K kms“^) 
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Figure 1. Schematic change in cold gas mass, H 2 mass, cold gas fraction (/gas = ^gas/(^gas + M*)) and molecular hydrogen fractions 
(/h2 = ^H 2/^gas) produced by our model when varying the stellar mass (left column), SFR (middle column) and scale radius (right 
column) of a galaxy. For each column the other two of three input galaxy properties of our model (stellar mass, SFR, and scale radius) 
are fixed. The fixed values are quoted in the upper panels. 


2011). One of the samples corresponds to the Great Obser¬ 


vatories Origins Deep Survey South (GOODS-S; Giavalisco 
et al.||20Q4 ), and the other to the UKIDSS Ultra-Deep Sur¬ 
vey (UDS) field (Lawrence et al. 20Q7| ). We used the pub¬ 
licly available GANDELS photometric catalogues compiled 
by Guo et al. (2013) and Galametz et al. (12013 ), and the 


redshift compilation by Dahlen et al. (2013|), with spectro¬ 


scopic redshift updates taken from Fadda et al. (2010) and 


Barro et al. (2014), to select our galaxies. We selected disc 


galaxies based on their Sersic index and excluded galaxies 
that contain an active galactic nucleus (AGN). We did not 
apply any selection based on stellar mass or SFR. The ex¬ 
act selection criteria and our methods to calculate stellar 
masses, SFRs, and sizes are described below. 

As our method requires us to know the scale radius of 
each galaxy, we extracted from these catalogues all those 
sources with good-quality morphological fits in the H band, 


as determined by van der Wei et al. (2012). We computed the 
circularized, effective radius Pe of a galaxy as Pe — Uca/ (b/a), 
where Ue is the effective radius along the major axis and 


b/a the ratio between the minor and major axis sizes and 
converted these to the scale length r* for a Sersic profile with 
Sersic index n. We converted all scale radii to a common rest- 
frame wavelength of 5000 A following 


van der Wei et al. 


(2014). Since our model is designed for disc galaxies, we 
discarded all galaxies with a Sersic index n > 2.5. 

We excluded galaxies known to contain an AGN from 
our sample, to avoid any bias in the derived star formation 
rates and stellar masses. To identify AGN, we used the 4Ms 


Ghandra X-ray catalogue for the GOODS-S field (Xue et al. 


2011 ), and an infrared power-law spectral energy distribu¬ 
tion analysis ( Gaputi|2013 ) in the UDS field, as the existing 
X-ray data is shallow ( Ueda et al.|2008 ). These AGN consti¬ 
tute only 1-2% of the GANDELS 0.5 < z < 3 samples. It is 
possible that other galaxies containing an AGN remained in 
our sample, but the AGN is likely weak as it is not identified 
in the X-ray data or with the IRAG power-law criterion. 

By applying these morphology and AGN cuts, we kept 
around 60% of the galaxies in the GOODS-S and UDS GAN¬ 
DELS catalogues. The distributions of the //-band magni- 
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Figure 2. if-band magnitude distributions of the original (white) 
and final samples (after the morphology cut; red) at redshifts 
0.5 < 2 ; < 3 for the GOODS-S (top) and UDS (bottom) fields. 


The Bruzual (2007) templates have a larger contribution 
from thermally pulsing asymptotic giant branch (TP-AGB) 
stars than the Bruzual & Chariot (2003) templates. Stellar 


masses calculated based on Bruzual 


(2007) are on average 
0.1 dex less massive than stellar masses calculated based on 
the Bruzual & Chariot (2003) models. Internal dust extinc¬ 
tion has been taken into account by convolving all Bruzual 
and Chariot templates with the Calzetti et al. (2000 | ) red¬ 
dening law, with Av values ranging from 0.0 to 3.0 (with a 
step of 0.1). 

SFRs are based on rest-frame UV fluxes corrected for 
extinction E(B-V) for each individual galaxy using the 


Calzetti et al. (2000) reddening law. UV fluxes were con¬ 
verted into SFRs following Kennicutt ( 1998a| for a Chabrier 
IMF. When IR photometry was available (for just a few per¬ 
cent of our sample) we calculated SFRs based on a combina¬ 
tion of non-extinction corrected UV and IR photometry. We 
derived rest-frame S/xm luminosities from the 24/im fluxes, 
and converted the 8/xm luminosities into total IR luminosi¬ 
ties following Bavouzet et al. (2008). SFRs were calculated 


following Kennicutt (1998a) and Bell et al. (2005). In com¬ 
parison with Popping et al. (2012) the CANDELS survey is 


closer to a mass-selected sample and (when IR photometry 
is available) is based on a more-reliable tracer of the SFR. 
As such, the results presented in this work supercede our 
previous efforts. 


3 RESULTS 

In this section we present galaxy SFRs and the inferred total 
cold gas, atomic hydrogen, and molecular hydrogen masses 
of galaxies in the CANDELS survey. We will focus on the 
gas masses, gas fractions, gas properties of galaxies and the 
evolution of the gas content of galaxies. 


tudes of the original and final samples at 0.5 < z < 3 for the 
two fields separately are compared in Fig.|^ From this Fig¬ 
ure, we can see that I) the GOODS-S sample reaches around 
a magnitude deeper than the UDS sample, as the former is 
based on both the deep and ultra-deep CANDELS data for 
the GOODS-S field; and 2) the percentage of galaxies lost 
because of the good-quality morphology cut and cut in Sersic 
index is not constant with iL-band magnitude. This cut still 
contains > 50% of the original galaxies for iL < 26, but the 
sampling becomes progressively more sparse at fainter mag¬ 
nitudes. The depth of the GOODS-S sample reflects itself 
in a relatively higher number of low-stellar mass galaxies. 
Nevertheless, our sample is likely to miss extended galaxies 
with low surface brightness. This becomes especially impor¬ 
tant for galaxies with low stellar masses at 2 ; > 2. For our 
model this means that we are biased against low surface 
brightness galaxies ( Guo et M]|20I3 ). 

We computed stellar masses using a multi-wavelength 
SED fitting to the CANDELS photometry (from the U 
to the 4.5 fim band), applying the Bruzual ( 2QQ7| ) templates 
fixed at the redshifts of the sample objects. We used a single 
stellar population and five exponentially-declining star for¬ 
mation histories (with r = O.I, 0.5, I.O, 2.0, and5.0 Gyr). 
For each star formation history we considered 24 possible 
templates, corresponding to ages between 0.05 and 5.0 Gyr. 


3.1 Galaxy SFR 

To fully appreciate the predictive power of our model it is 
crucial to place the inferred gas masses in their proper con¬ 
text. Our cold gas and H 2 estimates depend on the SFR of 
a galaxy to first order. In Figure we present the SFR of 
our galaxy sample as a function of stellar mass for differ¬ 
ent redshift bins. For clarity the relations are shown for the 
samples taken from the GOODS south field, the UDS field, 
and the combination of these fields. We find an increasing 
trend in SFR with stellar mass up to redshifts 2 : < 3. The 
median of this trend has been referred to as the ‘main se¬ 


quence’ of star-forming galaxies (e.g., Noeske et al. 2007 
Elbaz et al.||20II ) . The normalization of the trend between 


stellar mass and SFR decreases with time. At fixed stellar 
mass, galaxies in the redshift range 2.5 < 2 ; < 3.0 formed 
stars an order of magnitude more rapidly than in the red¬ 
shift range 0.5 < 2 ; < 1.5. At 0.5 < 2 ; < I.O the increasing 
trend of SFR with stellar mass has a slope of 0.66. Other 
studies have found a slope ranging from 0.6 to I at the same 
redshifts ( Pannella et al.|2009 Karim et al.|20Il Whitaker 


erar]|20I2[|20I4|). 

We did not apply any additional selection criteria be¬ 
yond those described in Section |2.2[ While many other 
studies only select galaxies on the ‘main sequence’ of star- 
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Figure 3. Galaxy SFR as a function of stellar mass in different redshift bins. The shaded regions show the log of the conditional 
probability distribution function P(SFR|M*), which represents the probability distribution of SFRs for fixed stellar masses, and the blue 
solid and dotted lines show the median fit and 2a deviation. The left column shows results for the GOODS-S sample, the middle column 
for the UDS sample, and the right column for the combined GOODS-S and UDS sample. The black line marks the mean trend of the UV 
+ IR based SFRs presented in |Barro et al.| ( [2014| > in the GOODS-S and UDS fields. The black dashed line mark the double-power-law 
fit to the relation between stellar mass and SFR presented in [Whitaker et ar] ( |2014| . 


formation, our sample also includes a fraction of galaxies 
with low SFRs. 


We compare our SFRs with SFRs from [Barro et al.j 
(2014 black line in Figure[^ and the double power-law fit to 
the stellar mass - SFR relation presented in [Whitaker et ah] 
(2014 black dashed line). We find that our mean trends for 
the stellar mass - SFR relation are in reasonable agreement 
with trends found in the literature. We note that IBarro et al.l 
(2014) and Whitaker et al. (2014) applied criteria to select 
actively star-forming galaxies, attempting to exclude more 
quiescent objects, whereas our selection criteria do not rule 
out quiescent objects. 


l.OM 0 yr ^). These galaxies with low SFRs are less ac¬ 
tively forming stars than counterparts at fixed stellar mass 


and make up a group of quiescent galaxies (Brennan et al. 
[2015[ ). For these stellar masses our sample is significantly 
incomplete (completeness < 50 per cent) below SFRs of 
log(SFR/Mo yr“^) = —0.5 at z = 2. The SFR limit to ob¬ 
serve galaxies increases at higher redshifts. This implies that 
at redshifts z > 2.0 the contribution by galaxies with low 
SFRs may be larger. 


3.2 Cold gas, HI, and H 2 content 


There is a large population of galaxies with stel¬ 
lar masses larger than ~ 10^° Mq and low SFRs (< 


We present the derived total cold gas masses of our galaxy 
sample as a function of stellar mass in different redshift bins 
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Figure 4. Total cold gas mass of galaxies as a function of their stellar mass in different redshift bins. The shaded regions show the 
log of the conditional probability distribution function P(Mgas|M*), and the blue solid and dotted lines show the median fit and one a 
deviation. Columns are as described in the caption of Figure 


in Figure We find that on average the cold gas mass of 
a galaxy increases with its stellar mass at all observed red- 
shifts with a sub-linear slope. Despite the clear increasing 
trend, we do find a significant amount of scatter. Note espe¬ 
cially a group of gas-poor galaxies with high stellar masses. 
This group of galaxies seems to have decoupled itself from 
the general increasing trend of gas mass with stellar mass, 
indicating that some physical process has either removed 
some of the cold gas or prohibits the cooling and/or ac¬ 
creting of gas onto the galaxy. These galaxies are likely on 
their way towards the red sequence. There is a shift in the 
normalization in the relation between gas- and stellar-mass 
with redshift. On average, high-redshift galaxies (especially 
at z > 2.0) are more gas-rich than galaxies with similar stel¬ 
lar mass at lower redshifts. At 0.5 < z < 1.0 the increasing 
trend of cold gas mass with stellar mass has a slope of 0.58. 


We present the derived H 2 masses of our galaxy sample 
in five redshift bins in Figure Similarly to the total cold 
gas mass, we find an increase in H 2 mass with stellar mass 
and a decrease in H 2 mass with decreasing redshift. We see 
a group of galaxies with high stellar masses and low H 2 
masses at < 2.5 in good agreement with a similar group 
of galaxies seen in the relation between cold gas and stellar 
mass (Fig. [^. The relation between H 2 mass and stellar 
mass is well defined and has a clear upper envelope over 
the entire range of stellar masses probed, with only minor 
stochastic appearance of galaxies above the envelope. The 
trend between H 2 mass and stellar mass is slightly steeper 
than for the total cold gas mass. At 0.5 < z < 1.0 the 
increasing trend of H 2 mass with stellar mass has a slope of 
0.65. Although the relations between stellar mass and cold 
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Figure 5. H 2 mass of galaxies as a function of their stellar mass different in redshift bins. The shaded regions show the log of the 
conditional probability distribution function P(Mh 2|^*), and the blue solid and dotted lines show the median fit and one a deviation. 
Columns are as described in the caption of Figure jsJPurple pentagons, circles, squares, and diamonds are literature values from [Paddij 
|et al.| ( [?010| >, [Tacconi et ^ ( |2010| >, [Tacconi et al.| ( [^13| >, and |Saintonge et al.| ( |2Q13| , respectively. 


and H 2 mass are very similar, there is no simple one-to-one 
mapping between cold gas mass and H 2 mass. 

As a comparison we have included direct measures of 
the H 2 mass of galaxies (through their CO luminosity) from 
the literature jTacconi et al.| 201^ Saintonge et al.| 2013 


Tacconi et al.||2013| in Figure pp 


The mean trend for the 
inferred H 2 masses lies below the observed sample. Our sam¬ 
ple contains 24000 galaxies covering four decades in stellar 
mass. The observational data points correspond to only a 
small number of galaxies in a very localized region of param¬ 


eter space. Our sample was not selected to purely consist of 
‘main sequence’ galaxies, driving the mean trends in inferred 
gas mass down with respect to a galaxy sample consisting 
of only ‘main sequence’ galaxies only. A close look at Figure 
shows that most of the directly-observed H 2 masses are 
in good agreement with the most-H 2 -massive (i.e., actively 
star-forming) galaxies at stellar masses M* > 

We hnd that the increasing relation between stellar mass 
and H 2 mass suggested by the direct observations continu¬ 
ous at least down to stellar masses of 10® M©. 


^ [Saintonge et al.| ( |2013| calculates the CO-to-H 2 conversion fac¬ 
tor as a function of metallicity using the prescription presented 


Genzel et al. 


( 2012 ). 


^ [Tacconi et al.| ( |2Q13| selected galaxies with a stellar masses ^ 
2.5x10^^ M© and star-formation rate^ 30 M©yr“^. This matches 
the most-actively star-forming galaxies in our sample. 









































































Gas evolution in CANDELS galaxies 9 



log M* [Mq] log M* [Mq] log M* [Mq] 


Figure 6. The HI mass of galaxies as a function of their stellar mass in different redshift bins. The shaded regions show the log of the 
conditional probability distribution function P(Mhi|M*), and the blue solid and dotted lines show the median fit and one a deviation. 
Columns are as described in the caption of Figure 


We show the HI mass of galaxies as a function of their 
stellar mass in different redshift bins in Figure Like for the 
cold gas and H 2 mass, the Hi content of galaxies increases 
as a function of stellar mass. The scatter around the mean 
trend is very wide, much wider than trends with cold gas 
mass and H 2 mass. At 0.5 < z < 1.0 the trend between HI 
mass and stellar mass has a slope of 0.54. 

The relations between stellar mass and cold gas, HI, 
and H 2 mass tightly follow the trend between SF and stellar 
mass (Figure]^. The well defined upper envelope of this 
relation for non-starburst galaxies is in good agreement with 
the upper envelope of the relation between stellar mass and 
H 2 mass. It is important to realize that in our model the 
H 2 mass is to first order a different representation of the 
SFR of a galaxy (through Asf), hence the similarity between 
the slopes of the relation of SFR and H 2 mass with stellar 
mass (0.66 for SFR versus 0.65 for the H 2 mass at z = 


0.5 — 1.0). The total cold gas mass and Hi mass are driven 
through a more complex combination of stellar mass, SFR 
and disc size, and therefore do not necessarily have a simple 
relationship with SFR. This is reflected in the difference in 
the slopes of the relations between stellar mass and SFR 
(slope of 0.66), cold gas mass (slope of 0.58), and HI (slope 
of 0.54). The inferred cold gas mass has a slightly shallower 
increase with stellar mass than the inferred H 2 mass and 
SFR. 


3.3 Gas fractions 

In Figure we show galaxy gas fractions as a function of 
stellar mass and redshift. There is a strong anti-correlation 
between galaxy stellar mass and gas fraction at z < 2.5. 
The cold gas fraction remains relatively constant at stellar 
masses M* ^ 10®'^ M© and drops rapidly at higher stellar 



























































10 G. Popping et al. 


. - GOODS-S + UDS 

1 .U I M ^-r 




1.0 


z = 2.5-3.0 


a 0.5 ^ 

44-1 



0 . 0 . 


8 9 10 11 12 

logM*[M0] 


- - 0.8 


- - 1.0 


- 1.2 


-1.4 


- 1.6 


- 1 .^ 


- 2.0 


- 2.2 



- 0.8 


- - 1.0 


- 1.2 


-1.4 


- 1.6 


- 1 . 


- 2.0 


- 2.2 


Figure 7. Cold gas fraction of galaxies (/gas = MH2+m+^ ^ 
as a function of stellar mass and redshift. The shaded regions 
show the log of the conditional probability distribution function 
F(/gas|-^*), and the blue solid and dotted lines show the median 
fit and one cr deviation. The black solid and dashed lines show 
the mean fit and one cr deviation to the predictions of [Popping, | 
[Somerville fc Tragerj ( [2QT4| . 


masses for galaxies in our highest redshift bin (2.5 < z < 
3.0). In the lower redshift bins the trend between stellar 
mass and cold gas fractions is roughly linear, mainly driven 
by a decrease in the cold gas fractions of galaxies with stellar 
masses Mq. The characteristic mass above which the 

gas fractions rapidly drop (M* ]^q9 .5 -io ]y[Q^ suggests that 

some physical process prevents the buildup of large gas reser¬ 
voirs in galaxies with a stellar mass above this characteristic 


mass, and is similar to the quenching mass scale (Kauffmann 
et al.|2Q03 ). We find the strongest evolution in galaxy cold- 


Figure 8. Relative molecular content as a function of stellar mass 
for different redshift bins. The shaded region shows the log of the 
conditional probability distribution function 1-^*); 

and the blue solid and dotted lines show the median fit and one cr 
deviation. The black solid and dashed lines show the mean fit and 
1 cr deviation to the predictions of [Popping, Somerville fc Trager[ 


erature values from [Daddi et al.[ ([2010J, 

Tacconi et al.[ ([2010|) 

[Tacconi et al.[ ([2013J, and[Saintonge et al. 

([2013J, respectively. 


gas fraction in intermediate mass galaxies (lO^^^^M©). At 
a fixed stellar mass, a galaxys cold-gas fraction decreases 
with time. 

We present /h 2 * = m^^+m ^ function of galaxy 

mass in Figurej^ This quantity is often used as an observable 
tracer of the gas content of galaxies. Similar to the total cold- 
gas fraction, there is an anti-correlation between H 2 and 
galaxy stellar mass, as well as redshift. High stellar mass 
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Figure 9. Relative atomic hydrogen content as a function 
of stellar mass for different redshift bins. The shaded region 
shows the log of the conditional probability distribution func¬ 
tion blue solid and dotted lines show 

the median fit and one a deviation. The black solid and dashed 
lines show the mean fit and 1 a deviation to the predictions of 
[Popping, Somerville &z, Tragerj ( |2014| ). 


galaxies have on average the lowest relative H 2 content, and 
/h 2 * decreases with redshift at a fixed stellar mass. The 
trend in fu 2 * with stellar mass is similar to that for total 
cold-gas fraction. 

We compare the results of our model with direct obser¬ 
vations (through CO) of /h 2 * in star-forming galaxies from 
Tacconi et al. ( |2Q10 ), Saintonge et al. (2013), and Tacconi 


et al. ( |2013| . Although direct observations of /h 2 * do not lie 


on top on the mean trend, they are in good agreement with 
the most massive and H 2 -rich galaxies in our sample. As dis¬ 


Figure 10. Molecular fraction (fH2 = MH 2 /Mgas) of the cold gas 
as a function of stellar mass for different redshift bins. The shaded 
region shows the log of the conditional probability distribution 
function and the blue solid and dotted lines 

show the median fit and one a deviation. 


cussed when we described Figure the difference between 
the mean trend and observations is driven by different selec¬ 
tion criteria. Tacconi et al. ( 2013| developed a formalism to 
correct gas fractions for this selection effect and found gas 
fractions of ~30 per cent at = 1 — 1.5 for galaxies with 
stellar masses ~ 10^^ Mq, in reasonable agreement with our 
inferred gas fractions. 

We present the relative amount of atomic hydrogen in 
galaxies (/hi* = ^ function of stellar mass for 

different redshift bins in Figurej^ Similar to the previous fig¬ 
ures, /hi* decreases with stellar mass. Below stellar masses 
of ~ 10® M 0 the Hi mass of galaxies exceeds the stellar 
































































12 G. Popping et al. 



Figure 11. Evolution of the mean trends in /gas, /h 2 *, /hi*, and 
/h 2 as a function of stellar mass. Note that the mean trend in 
/hi* with stellar mass remains constant with time, whereas the 
other three properties do evolve. 


mass, whereas at higher stellar mass the HI mass is less 
than the stellar mass. Interestingly, there does not seem to 
be a strong evolution in the trend between /hi* and stellar 
mass with time. A similar transition mass from Hi to stellar 


mass dominated is seen for local galaxies (e.g., |Cortese et ah 
Tluang et ^|2012). 


2011 


In Figure |10| we present the molecular fraction of the 
cold gas (fH 2 = MH 2 /Mgas) as a function of stellar mass 
for different redshift bins. This is a good measure of the 
fraction of gas available for SF as a function of time. The 
relation between fH 2 and stellar mass is relatively flat at 
z > 2.5. We find an increasing trend between fH 2 increasing 
with stellar mass towards lower redshifts. We find that fH 2 
increases with stellar mass up to a stellar mass of ~ 10^° M© 
and decreases at higher stellar masses in our lowest redshift 
bins (0.5 < z < 1.0). The average H 2 fraction of galaxies 
decreases with time at fixed stellar mass (only at M* > 
10^° M© do galaxies at 2.0 < z < 2.5 have higher fH 2 than at 
z > 2.5). This result clearly shows that both the cold gas and 
H 2 reservoirs decrease with time, although not necessarily 
at the same rate. It is important to note that the relation 
between stellar mass and H 2 fraction has a large scatter, 
especially at low stellar masses and lower redshifts. There is 
no well defined drop in H 2 fraction at a characteristic stellar 
mass and/or time. 

We show the mean trends of the previous four Figures 
again in Figure pT] with the aim of focusing on the evolution 
of these trends. We find the remarkable result that /gas, 
/h 2 *, and /h 2 all decrease as a function of time, whereas 
/hi* remains nearly constant. Although the gas content of 
galaxies at fixed stellar mass decreases as a function of time, 
the changed partitioning of cold gas in HI and H 2 results in 
a constant fraction of HI in galaxies with fixed stellar mass. 
We discuss this further in Section [S] 

The CANDELS results suggest a gradual evolution in 
the cold gas fraction and the relative molecular hydrogen 
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Figure 12. Galaxy SFR as a function of stellar mass for different 
redshift bins. The colour map gives the cold gas fraction of the 
galaxies. Each galaxy is counted only once to create the hexagons. 


content of galaxies at 0.5 < z < 3, independent of stellar 
mass. 


3.3.1 Model comparison 

As a comparison we included predictions from the semi- 
analytic model presented in Popping, Somerville &: Trager| 
(2014) in Figures and|9 We ran the models in the ap¬ 
propriate redshift bins and applied a CANDELS selection 
criteria to the model output {Hab < 26). We only selected 
disc-dominated galaxies (bulge-to-total ratio less than 0.4). 

Although the results presented in this paper are merely 
model predictions, they are based on actual observables 
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Figure 13. Galaxy SFR as a function of stellar mass for different 
redshift bins. The colour map gives the molecular fraction of the 
cold gas. 


rather than introducing uncertain recipes describing the 
physics acting on galaxies. A comparison between the re¬ 
sults presented in this work and semi-analytic predictions 
can shed light on the successes and failures of theoretical 
models. We acknowledge that the results presented are ef¬ 
fectively a different representation of the phase space de¬ 
scribed by stellar mass, SFR, and galaxy size. One could 
therefore also compare the observed values for these quan¬ 
tities to semi-analytic model predictions. Nevertheless, fo¬ 
cusing on the gas masses of galaxies and potential discrep¬ 
ancies between the results presented in this work and semi- 
analytic model predictions may provide a more direct hint 
of the baryonic physics that need to be changed and/or in¬ 


cluded for models to reach better agreement with observa¬ 
tions. Somerville & Dave (2014) showed that semi-analytic 
models and hydrodynamic models predict too little SF in 
galaxies at redshifts z — 1 and z — 2 (for a potential solu¬ 
tion see |Henriques et al.|2014|[White, Somerville &: Ferguson| 


2015). Our method has the potential to assess whether this 


is due to too low gas masses in general or too little H 2 . 

The semi-analytic model predictions are in good quali¬ 
tative agreement with our inferred cold gas fractions at red- 
shifts z = 0.5toz = 2.5 (Figure [^. We predict up to 10 per 
cent more cold gas in galaxies with stellar masses larger than 
10 ^° M 0 than the semi-analytic model. At higher redshift 
the semi-analytic model predicts gas fractions higher than 
found in this work for galaxies with stellar masses larger 
than 10 ^° M 0 . 

Semi-analytic model predictions are in good qualita¬ 
tive agreement with our inferred H 2 masses at redshifts 
1 < z < 1.5 (Figure]^. At lower redshifts the semi-analytic 
model predicts much less molecular hydrogen than the re¬ 
sults presented in this paper, whereas the model predict on 
average 1.8 times less molecular hydrogen at higher redshifts 
(/h 2 * is approximately 10 per cent less). In the highest- 
redshift bin the model overpredicts /h 2 * compared to the 
inferred fractions at M* > 10 ^° M 0 , whereas it underpre¬ 
dicts the relative amount of H 2 in galaxies at lower stellar 
masses. 

The Popping et al. semi-analytic model on average 
predicts higher values for /hi* at stellar masses less than 
10^° M 0 than the results presented in this work. The model 
also finds that at redshifts z = 3 — 1 the amount of Hi 
in galaxies at fixed stellar mass remains relatively constant. 
We will further discuss the implications of this comparison 
in the discussion. 


3.4 Gas properties on the stellar mass — SFR 
diagram 

In Figure we show the cold gas fraction of galaxies on 
the stellar mass-SFR relation. We find that to first order, 
galaxy cold gas fractions are well parameterized by their 
stellar mass and SFR. It is important to realize that both 
the SFR and the stellar mass of a galaxy are key ingredi¬ 
ents in our model that set a galaxy’s cold gas mass and its 
partitioning into atomic and molecular hydrogen. Galaxy 
gas fractions stay constant parallel to the relation between 
SFR and stellar mass, increasing towards higher SFR and 
decreasing with stellar mass. The upper envelope in the stel¬ 
lar mass-SFR diagram is populated by the most gas-rich 
objects at a given stellar mass. Although counterintuitive, 
galaxies with the highest SFRs do not necessarily have the 
highest gas fractions. The highest gas fractions are found 
in galaxies with low stellar masses and high SFRs. Galax¬ 
ies with high stellar mass and low SFRs have the lowest 
gas fractions. At a fixed position in the stellar mass-SFR 
diagram galaxy cold gas fractions decrease with time. 

We explore the H 2 fraction of galaxies in the stellar 
mass-SFR plane in Figure Molecular gas fractions in¬ 
crease with increasing SFR. The most actively star-forming 
galaxies have highest molecular gas fractions, and the least- 
act ively star-forming galaxies have the lowest fractions. 
Molecular fractions cannot be parameterized by a simple 
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combination of stellar mass and SFR as the total cold-gas 
fractions. 

More information can be obtained by comparing Fig¬ 
ures!^ and These figures clearly show that large cold- 
gas reservoirs do not necessarily lead to large molecular-gas 
reservoirs and active SF. This is especially evident for galax¬ 
ies with low stellar mass (~ 10® M©) and SFR ~ 1 M© yr“^ 
at z < 1.5. These galaxies are dominated by atomic gas. De¬ 
spite being rich in cold gas, only approximately 25 per cent 
of their cold gas is molecular and available to form stars. 
These objects were below the observation limit of the galax¬ 
ies in Popping et al. (2012). The deep near-IR photometry of 


the CANDELS survey allows us to study these atomic-gas- 
dominated low-mass galaxies. Observations of local galaxies 
show a similar decrease in H 2 fraction with decreasing stel¬ 
lar mass (Saintonge et al. |2011|[Boselli et al.|2014[ |BothwelT| 


et al.|[Ml4 ), with on average even lower H 2 fractions than 

the galaxies in our sample. 

The most actively star forming galaxies (~ 10^^ M© 
and SFR ~ 100M©yr“^) al z < 2.5 have relatively small 
cold gas reservoirs (fgas < 0.5), but most of this cold gas is 
molecular. Galaxies with high stellar mass and low SFRs are 
nearly empty of cold gas, but the molecular fraction of this 
gas is approximately 0.5 and can account for some residual 
star-format ion. 


3.5 Gas properties on the stellar mass — size 
relation 

We present the cold-gas and molecular fraction of galax¬ 
ies on the stellar mass-size relation in Figures and |15| 
Within our model the scale radius of a galaxy sets the sur¬ 
face densities that control the molecular fraction and the 
total mass of cold gas in that galaxy. At fixed stellar mass 
the cold gas fractions of star-forming galaxies increase with 
increasing size. The molecular fraction of the cold gas, on 
the other hand, decreases with increasing size. For a fixed 
scale radius the cold-gas fraction and molecular fraction of 
a galaxy decrease with time. 

The largest variations in molecular fractions are found 
in low-stellar-mass galaxies (M* < 10^° M©). At higher stel¬ 
lar masses we find differences in molecular fraction of ap¬ 
proximately 25 per cent. Although less dramatic than for 
lower-mass galaxies, such variations are still significant. 


3.6 Star-formation efficiencies 

The star-format ion efficiency (the ratio between SFR and 
total cold gas mass (Hi + H 2 ), i.e., SFE = SER/Mgas) is a 
good measure of how efficiently galaxies turn cold gas into 
stars|^We show the SEE of galaxies on the stellar mass-SER 
and stellar mass-size diagrams in Eigures and pT| 

We find that galaxy SEEs increase with increasing SER 
and decrease with time. The SEE of galaxies that populate 


^ We note that the star-formation efficiency of galaxies is often 
defined as the ratio between SFR and molecular hydrogen mass 
H 2 (i.e., SFE = SFR/Mh 2 )- Our adopted definition of the ratio 
between SFR and total cold gas mass is slightly different and 
important for simulations that do not discriminate between HI 
and H 2 . 
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Figure 14. Galaxy scale radius as a function of stellar mass for 
different redshift bins. The colour map gives the cold gas fraction. 


the upper envelope in the stellar mass-SFR relation changes 
as a function of stellar mass. At high stellar masses, galax¬ 
ies are more than twice as efficient in forming stars out of 
their total cold gas reservoirs (Hi + H 2 ) than at low stellar 
masses. This shows that the position of a galaxy on the stel¬ 
lar mass-SFR plane is being driven by the amount of cold 
gas available and by the ability of the gas to form molecules 


and stars ( 

Saintonge et al. 2012 Genzel et al. 2014 Sar- 

gent et al. 

2014). The group of galaxies with low SERs at 


high stellar masses is more than twice as inefficient in form¬ 
ing stars than their counterparts with the same stellar mass 
that are actively forming stars. Overall, the SFE of galaxies 
decreases with time. 

The SEE of galaxies decreases with increasing scale ra- 
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Figure 15. Galaxy scale radius as a function of stellar mass for 
different redshift bins. The colour map gives the molecular frac¬ 
tion of the cold gas. At fixed stellar mass more compact galaxies 
have higher molecular fractions. 


Figure 16. Galaxy SFR as a function of stellar mass for different 
redshift bins. The colour map gives the star-formation efficiency 
(SFE = SFR/Mgas) of the galaxies. 


dius. This is especially clear for galaxies with low stellar 
mass, but is also apparent for the most massive galaxies in 
our sample. This result is in good agreement with our find¬ 
ings in Figure p~5l Compact galaxies have higher molecular 
fractions and are therefore more efficient at forming stars. 


3.7 Global evolution of cold gas content of 
galaxies 

The evolution of /h 2 * with redshift has been studied in sev¬ 
eral different CO surveys and through dust continuum esti¬ 
mates (Geach et al. ] |2011[ |Tacconi et al.|[MTQ Daddi et al, 


2010 

IMagdis et al.|2012||Tacconi et al.|2013| 

Saintonge et al. 

2013 

IBauermeister et al.||2013| |Genzel et a 

.||2014| IScoville 


et al. |2014 ). These surveys suggest /h 2 * increases with in¬ 

creasing redshift up to z ~ 2 and flattens or even declines 
at higher redshifts ( Saintonge et al.|2013 ). 

We present the evolution in fu 2 * for galaxies with dif¬ 
ferent selection criteria in stellar mass and SFR in Fig. 
(top row). Our method allows us to infer the gas content for 
a large number of galaxies covering a wide range in physical 
properties (stellar masses, SFRs). This makes it ideal to in¬ 
vestigate the effects that different selection criteria have on 
cold gas evolution trends inferred from direct observations 
of the CO emission line and the dust continuum of galaxies. 
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Figure 17. Galaxy scale radius as a function of stellar mass for 
different redshift bins. The colour map gives the star-formation 
efficiency (SFE = SFR/Mgas) of the galaxies. 


We find a decrease in the relative H 2 content of galaxies from 
z ~ 3 to z ~ 1. There is a minor increase in the relative H 2 
content (and other fractions) from 2 := 1 . 0 toz = 0.5 when 
selecting galaxies with stellar masses larger than 10^ M©. 
Above z = 2.0 the evolution of /h 2 * becomes less steep. This 
observed evolution is similar to the results presented in |Sain-| 


tonge et al. (2013). The normalization and the shape of the 


evolution in /h 2 * depend on the selection criteria applied. 
We find the best match with the current direct observations 
when selecting only galaxies with SFR > 30 Mq yr ^ for 
galaxies with stellar masses log (M* /M©) > 10.0, similar to 
the selection limits of the observed samples. When loosen¬ 
ing the selection criteria for SFR the normalization of /h 2 * 


decreases. These galaxies with low SFRs are currently typ¬ 
ically not accounted for in direct surveys of the molecular 
content of galaxies through either CO or FIR sub-mm con¬ 
tinuum studies. When including galaxies with lower stellar 
masses we find that not only the normalization in the evo¬ 
lution of /h 2 * increases, the slope of the decline in /h 2 * at 
z ^2 also becomes stronger. We find no changes in the evo¬ 
lution of /h 2 * when selecting galaxies with stellar masses 
more massive than 10^ solar masses and changing the SFR 
criteria. 

The relative amount of atomic hydrogen in galaxies 
slightly decreases with time in galaxies with stellar masses 
larger than 10® M© and remains constant in galaxies with 
stellar masses more massive than 10^® M©. 

The mean molecular fraction fH 2 of the cold gas in 
galaxies gradually decreases over the entire redshift range 
probed. This is in good agreement with the results presented 
in Figure where we showed that the mean molecular 
fractions as a function of stellar mass decreases by approxi¬ 
mately 25 per cent for the lowest mass galaxies from redshift 
3 to 0. 


4 MODEL ASSUMPTIONS AND 

APPLICABILITY TO THE CANDELS 
SAMPLE 


The presented approach is based on several assumptions that 
may not always hold. In this section we discuss the main 
assumptions and how they affect our inferred gas masses. 

We converted all the sizes to a common restframe wave¬ 
length of 5000 A. Originally they were based on H band 
photometry. Across the 2 : = 0.5-3 redshift range the H-hand 
photometry spans a restframe size from the B band to the 
I band. Inferred gas masses based on the 5000 A restframe 
morphology are slightly higher than based on the original 
iF-band morphology. When adopting the 5000 A restframe 
morphology, less than 5% of the galaxy sample is more than 
10 % more gas rich than when adopting the iF-band mor¬ 
phology. Only a handful of objects is more than 20% more 
gas rich when adopting a 5000 A morphology, rather than 
the iF-band morphology. 

A key assumption in our method is an exponential dis¬ 
tribution of stars and gas in the galaxy discs. SF takes place 
in molecular clouds (local clumps in the disc not follow¬ 
ing an exponential distribution) which could lead to a lo¬ 
cal underestimation of the cold gas surface density. In this 
work we are averaging over ~kpc scales and at our redshifts 
of interest and at these scales the exponential disc model 
is a good approximation to the stellar distribution of star¬ 
forming galaxies on galactic scales (Wuyts et al. |2011| |. Fur- 
thermore, analyses of local galaxies have revealed an expo¬ 
nential distribution of the cold gas in star-forming galaxies 


Leroy et al.|2008l Bigiel fc Blitz|20l2 Kravtsov|2013 Wang 


et al.||2014 ). 

Our model assumes a constant ratio between the gas 
and stellar disc scale lengths (xgas = 1-7), based on fits to 
the galaxy disc profiles in Leroy et al. (2008). This number is 


in good agreement with typical ratios of 1.5-2.0 between Hi 
and optical disk sizes ( Verheijen Sancisi||2001 ). Kravtsov 
(2013) finds a Xgas = 2.6, when normalizing gas profiles to 
a radius Pn obtained using abundance matching arguments. 
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Figure 18. Redshift evolution of the relative H 2 content 
row), gas fraction (second row), relative HI con¬ 
tent (third row), and molecular hydrogen fraction (bottom row) 
for galaxies with stellar masses in the range M* > 10^ Mq (left 
panels), M* > 10^° Mq (right panels), and SFR larger than 0 
(red), 5 (green), 10 (orange), and 30 (blue) MQyr“^. Different 
markers represent the following datasets: upwards pointing tri¬ 
angles: galaxies from the THINGS survey ( [Leroy et al.||2008|> at 
2 : = 0.0, pentagons: galaxies from the COLD GASS survey (jSain- 
tonge et al. 2011| > at 2 ; = 0.0, squares: galaxie s fromjBauermeister 
et al.j ( [2013 1 at 2 : ~ 0.4, stars: galaxies from jGeach et al.|(|2011| 
at 2 ; ~ 0.4, eireles: galaxies from the PHIBBS survey (jTacconij 
jet al.|[M^ at 2 : ~ 1.2 and 2 ; ~ 2.2, downwards pointing trian¬ 
gles: galaxies at z ^ 1.2 and 2 ; ~ 2. 2 fromjTacconi et al. ( |2010| >, 
diamonds: galaxies at 2 : ~ 2.2 from ( jSaintonge et al.|2013 ). Error 
bars show the 80 percent confidence interval for each observa¬ 
tional sample at the given redshift bin. 


Furthermore, the value for Xgas may vary as a function of 
redshift or galaxy properties. Varying Xgas results in negli¬ 
gible differences in inferred cold gas and H 2 mass as long as 
Xgas > 1- Decreasing Xgas to values less than one leads to 
more significant differences in the inferred gas masses. When 
adopting the extreme case of Xgas = 0.5 we find that the in¬ 
ferred gas and molecular masses are lowered by 0.25 and 0.1 
dex, respectively. Observations of the sizes of the CO discs 
at z ~ 1.2 and z ^ 2.2 (supposedly tracing the molecular hy¬ 
drogen) do not support a scale length ratio of Xgas <1.0 (e.g. 
Tacconi et al.|2010 2013). Berry et al. (2014) found that ir 


semi-analytic models HI discs have to be more extended 
than the stars in order to reproduce Damped Lyman-alpha 
properties ( Schaye|2001 ). 

A key ingredient in our model is the partitioning of cold 
gas into a molecular and atomic component using a pressure- 
based recipe presented in Blitz & Rosolowsky (2006). This 
recipe was calibrated based on observations in local galax¬ 


ies (Blitz & Rosolowsky 2006 Leroy et al. 2008) and has 


not yet been constrained at our redshifts of interest. An 
increased UV background or lower gas metallicities could 
suppress the formation of molecules on dust grains and the 
ability of the molecular hydrogen to self-shield. Under these 
conditions our model would require a larger inferred cold gas 
reservoir to support the observed SFRs. We anticipate that 
future high-resolution observations with ALMA and SKA 
will be able to assess the validity of these locally calibrated 
relations at high redshifts. An alternative approach would 
be to use a molecular hydrogen prescription that is also a 


function of gas phase metallicity (e.g.. 

[Krumholz, McKee 

& Tumlinson||2009 Gnedin & Kravtsov 

[ 2011 ). Similarly to 


the pressure based algorithm, these prescriptions have also 


only been tested against local observations (e.g. Fumagalli 
Krumholz V Hunt|2010 ). A significant downside of applying 
this approach is that gas-phase metallicities are only avail¬ 
able for a limited number of galaxies at our redshifts of in¬ 
terest. A different approach could be to use the fundamental 
metallicity relation ( Mannucci et al.]|2009 ) to estimate the 
gas-phase metallicity of galaxies as a function of their SFR 
and stellar mass. We opted not to use this approach, as this 
will introduce an additional empirical step and source of un¬ 
certainty. Moreover, the calibration of gas phase metallicity 
measurements is highly uncertain. 

We assume an adapted version of the SF relation pre¬ 


sented in Bigiel et al. (2008) to allow for an increased star 


formation efhciency in high-surface-density environments. 
When not accounting for the increased star-formation efh- 
ciency (i.e. Vsf = 0 ), the presented method predicts larger 
values for the gas content of galaxies with SFR/(27rr*) 
lMQyr“^kpc“^ by roughly 10 percent and SFR/(27rr*) 

10 MQyr“^kpc“^ by roughly 75 percent (less than 1 per¬ 
cent of our sample has a SFR surface density larger than 
SFR/(27rr2) 10 Moyr-^kpc-^). 

To estimate the mid-plane pressure on the disc we adopt 


the prescription presented in Elmegreen (1989). We esti¬ 


mate the ratio between the gas and stellar vertical veloc¬ 


ity dispersion following Fu et al. (2010), based on empirical 


scalings for nearby disc galaxies. Observations suggest that 
star-forming galaxies at z ~ 2 are “puffy” with high velocity 
dispersions ( Forster Schreiber et al.|2009 ). Depending on the 
actual distributions of the giant molecular clouds and star¬ 
forming clumps, this may change the shielding properties of 
the molecular hydrogen. 

An additional source of uncertainty lies in the conver¬ 
sion from CO emission to H 2 mass (see [Bolatto, Wolfire 
Leroy 2013| for a review) or the conversion from sub-mm 
dust continuum to gas mass. Our method is calibrated on 
H 2 masses that rely on the aforementioned conversions. Our 
method therefore includes these uncertainties. 

In (e.g.. Popping et al.|2Ql^ , Figure 3) we showed that 
propagating the uncertainty in the input parameters (stel¬ 
lar mass, SFR, and size), the typical systematic uncertainty 
in inferred gas masses is approximately 0.3 dex. Many of 
the uncertainties cannot be well quantified. Combining all 
the concerns discussed above as well as the concerns dis¬ 


cussed in Popping et al. (2012 Section 3.3), and taking the 
errors in the input parameters into account, we estimate 
an uncertainty in the inferred gas masses for our method is 
approximately 0.7 dex. 
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5 DISCUSSION 


5.1 Galaxy gas fractions 

The gas fraction of a galaxy is set by the competition be¬ 
tween the inflow, outflow, and consumption of cold gas in 


the galaxy (Bouche et ah 2010 Dave, Finlator & Oppen- 


heimer 2011 

L Dekel & Mandelker 2014). The cosmic SER 

peaks at z ~ 2 (Hopkins & Beacom||2006 Madau & Dick- 

inson||20I4). 

, and analytic calculations and cosmological hy- 
al simulations predict large gas accretion rates 
)s at this epoch (Birnboim & Dekel||2003 Keres 

drodynamic 
onto galaxie 

et al.|2005 ' 

2009 Dekel & Birnboim|2006 Dekel et al.|2009). 


Our indirect gas measures provide an ideal sample to con¬ 
strain which of the aforementioned processes dominate this 
competition at different cosmic times. We find that gas frac¬ 
tions decrease most rapidly at z < 1.5, although a gradual 
decrease is present at all redshifts studied in this sample. 
These results suggest that the accretion of gas onto galax¬ 
ies cannot keep up with the consumption and/or outflow of 
gas, especially at redshifts z < 2 (see also [Popping et al. 


2012). These results are in good qualitative agreement with 


predictions from semi-analytic and semi-empirical models 


(e.g., 

Obreschkow & Rawlings 

|200^ Lagos et al.pOlll |Pu 

et al. 

|20I2| [Popping, Somervi 

le & Trager[[20I4| [Popping, 

Behroozi & Peeples||20I5|). Besides the relative gas con- 


tent, the molecular fraction of the cold gas also decreases 


with time, most rapidly at z < 2 (Saintonge et al. 2013). 


The characteristic density of the cold gas decreases ( Pop¬ 
ping et al.|2014 ), allowing less gas to self-shield and become 


molecular. This indicates that the physical process that sup¬ 
presses star-format ion with time is at least two-fold: galaxies 
run out of gas as well as molecules, but not necessarily at 
the same rate. For the most massive galaxies in our sample 
this two-fold process is already taking place at z ~ 3, well 
before the peak in cosmic SFR and before the accretion of 
cold gas onto galaxies slows due to the increasing dominance 
of dark energy driven expansion. This is an important con¬ 
straint for models of galaxy formation that include H 2 -based 
star-format ion recipes. 

The rapid decrease in cold gas fractions below z < 1.5 
is supported by direct observations of the CO emission in 
galaxies (Figure Tacconi et al. 2013 Saintonge et al. 


2013). These observations indicate that the relative H 2 con¬ 
tent of galaxies remains flat above redshifts of z ~ 2 and 
gradually decreases at lower redshifts. The slope and nor¬ 
malization of this trend is in part set by the selection crite¬ 
ria. Depending on the selection criteria in stellar mass and 
SFR chosen, the evolution of /h 2 * at redshifts z > 2 can 
change. 


Popping et al. (2012) found for a sample of galaxies 


taken from COSMOS at redshifts 0.5 < z < 2 that the cold 
gas fraction of galaxies with stellar masses of ~ 10^°'^ M© 
evolves from fgas ~ 1 at 1.75 < z < 2 to fgas ~ 0.2 at 
0.5 < z < 0.75. The authors found that less-massive galax¬ 
ies showed hardly any evolution in their gas fractions. This 
difference suggested that massive galaxies become gas-poor 
earlier and quicker than less massive galaxies; this behaviour 
was characterized as downsizing in gas content. Further¬ 
more, the authors found a characteristic stellar mass below 
which the gas fraction of galaxies rapidly drops. This char¬ 
acteristic mass decreased towards lower redshifts. We find a 
similar characteristic stellar mass at redshift 2.5 < z < 3, 


but the drop in the relation between cold gas mass and stel¬ 
lar mass disappears towards lower redshifts. Also, there is no 
strong difference in the rate at which the highest mass galax¬ 
ies run out of cold gas compared to the least massive galaxies 
in our sample (see FigurepT]). We remind the reader that the 
CANDELS survey allows us to study much fainter objects 
and is not biased towards the most star-forming objects. 
We ascribe the conclusion reached in Popping et al. (2012) 
to the selection of only the most actively star-forming, and 
therefore most gas rich, galaxies at stellar masses less than 
10^'^ M©. Our results suggests a gradual evolution in the 
cold gas fraction and the relative molecular hydrogen con¬ 
tent of galaxies at 0.5 < < 3, independent of stellar mass. 


star-forming objects 


Popping et al. (2012 


5.2 Extended parameter space and selection 
effects 

With the technique used in this work we can easily infer the 
cold gas masses for large samples of galaxies, covering a wide 
range in stellar mass, SFR, and redshift. As discussed above, 
this allows us to evaluate the impact of selection effects on 
smaller surveys of directly measured gas masses and to ex¬ 
tend observed relations to stellar mass, SFR, and redshift 
regimes where gas measure are hard to obtain directly. 

CO observations of actively star-forming galaxies have 
revealed a linear relation between H 2 mass and stellar mass 
(e.g. Saintonge et al.|2013 Tacconi et al.|2Q13 ). Our results 
suggest that this relation was already in place at z = 3, 
when the Universe was less than 3 Gyr old and extends 
down to stellar masses of at least 10® M©. A similar relation 
(although less well defined) is found for the total cold gas 
mass and HI mass. The trend between stellar mass and Hi 
mass has larger scatter than the trend between stellar mass 
and cold gas mass and H 2 mass. Even though the uncertain¬ 
ties in the gas masses are significant, these extended scaling 
relations suggest that the balance between the consumption 
and in- and outflow of gas in galaxies is already at place for 
galaxies with low stellar masses. 

In Eigure we explored the direct effect of changing 
selection criteria on the evolution of gas fractions in galaxies. 
We found no qualitative difference between different selec¬ 
tion criteria. The fraction of gas, HI, and H 2 in galaxies 
decreases gradually from 2 ; = 3 to 2 ; = 0.5. The absolute 
values of the gas fractions change by a factor of a few when 
changing the selection criteria. We find that for a sample of 
galaxies with stellar masses more massive than 10® M© the 
inferred cold gas. Hi, and H 2 fractions of galaxies remain 
nearly identical when galaxies with a reasonable amount of 
SE are selected (SERs of a few solar masses per year). Eu- 
ture surveys with a stellar mass limit of 10® M© should be 
able to robustly constrain the evolution of the gas content 
of galaxies contributing most to the stellar mass budget and 
stellar growth in our Universe. 

We found good agreement between direct observations 
of the H 2 mass of galaxies and our inferred H 2 masses 
for the galaxies with the highest SER in our sample (with 
SER > 20 — 30M©yr“^). Our results indicate that at the 
stellar masses probed by direct surveys, there is also a sub¬ 
stantial population of galaxies with H 2 masses less massive 
than 10^® M©. This group of galaxies should be accounted 
for when trying to make more generalized statements, espe¬ 
cially when focusing on galaxies with large stellar masses. 
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5.3 A constant relation between HI mass and 
stellar mass 

We found hardly any evolution in the relation between Hi 
mass and stellar mass, whereas the relation between cold gas 
mass and H 2 mass with stellar mass does evolve. [Popping 


Behroozi & Peeples (2015) found a similar result when cou¬ 


pling an abundance matching model with the semi-empirical 
approach to infer gas masses adopted in this work. Semi- 
analytic, semi-empirical, and hydrodynamic models have 
also found a constant relation between HI mass and stellar 
mass remains constant with time (for galaxies with stellar 
masses larger than 10^ M©), just as the HI mass function out 


to 2 : = 2 ([Dutton, van den Bosch & Dekel|2010[ Lagos et al. 

20111 [Eu et al.ll 

|2012| Dave et al.[[2013[ Popping, Somerville 

& Trager|2014[ 

Popping, Behroozi & Peeples|2015 ). 


Although already predicted by different models, this is 
the first time a model used actual observations as input to 
find a constant H I-to-stellar mass ratio. The relatively weak 
evolution in HI mass is driven by an apparent self-regulation 
that naturally arises in galaxies. As gas is consumed over 
time and galaxies grow as a function of stellar mass, the cold 
gas surface density decreases which naturally causes the HI 
fraction of the cold gas to increase. Although this process 
may prevent a rapid drop in the Hi reservoirs of galaxies, 
by no means does this process imply that the HI mass at 
a fixed stellar mass should remain constant. There is a fine 
balance between the amount of cold gas being consumed, 
heated and/or expelled, and the partitioning of the remain¬ 
ing gas into an atomic and molecular component that sets 
the constant Hi mass. 

We are aware that this result should be placed within 
the context of the model framework adopted in this (and 
other) work. If the adopted recipes for SF and the parti¬ 
tioning of cold gas are not valid for galaxies at z > 0 this 
could also dramatically change our conclusions. We therefore 
look forward to surveys with the SKA pathfinders probing 
the HI mass of galaxies out to at least z = 1. If indeed the 
H I-to-stellar mass fraction of galaxies remains constant, it 
should be feasible to observe the HI mass of M* galaxies at 
z = 1-1.5 ( Carilli Rawlings|2QQ4 ). 


5.4 Comparing predictions by a theoretical model 
to inferred gas masses 

Semi-analytic models of galaxy formation start from a cos¬ 
mological framework for structure formation, and include 
physical recipes to describe the baryonic physics acting on 
galaxies (such as the cooling and accretion of gas, star forma¬ 
tion, and outflows). In this work we started from observed 
quantities and converted them into other quantities using 
empirical scaling laws. Although many of the ingredients 
are the same for the two approaches, their different nature 
makes them very complementary. Any tension between pre¬ 
dictions from semi-analytic models and this work can shed 
light on the successes and failures of theoretical models. 


We find that the Popping, Somerville & Trager (2014) 
semi-analytic model predicts three times less molecular hy¬ 
drogen than our inferred gas masses suggest at redshifts 
0.5 < z < 1.0. The semi-analytic model predicts on av¬ 
erage about 2 times less molecular hydrogen at redshifts 
1.0 < z < 2.5. The same holds for galaxies at redshifts 


2.5 < z < 3 with stellar masses less than 10^° M©. Simi¬ 
larly, we find that the semi-analytic model predicts less cold 
gas in galaxies at redshifts 0.5 < z < 3 than our inferred gas 
masses suggest. 

Although the difference between semi-analytic model 
predictions and our inferred gas masses is not huge (and 
well within the uncertainty range of the method used in this 
work), it fits well within a broader picture where theoretical 
models predict too little SF for galaxies at redshifts 1 < 
z < 3 ( [Somerville Da'^|2014 ). If the molecular hydrogen 
masses of modeled galaxies are too low for their stellar mass, 
the SFRs of those galaxies will be too low as well when a 
molecular gas based star-format ion recipe is used. 


Popping, Behroozi & Peeples (2015) extended the sub¬ 


halo abundance matching approach with recipes to infer cold 
gas masses. They found that the same semi-analytic model 
as used in this work predicts less molecular hydrogen and 
cold gas at redshifts 1 < 2 : < 3 than their extended sub¬ 


halo abundance approach suggests (see also Somerville & 

Dave 2014 White, Somerville & Eerguson 

2015). A com- 

parison with different semi-analytic models ( 

Obreschkow & 


making this a more general problem. 


5.5 What drives the star-for mat ion efficiency of 
galaxies 

We have shown that, when using our model, the molecular 
fraction of the cold gas is a very important quantity that 
cannot be neglected. Despite large total cold gas masses, 
low molecular fractions can suppress the formation of stars 
and result in low SFE[^ Indeed the SFEs of galaxies on the 
relation between SER and stellar mass increases with stellar 
mass ( Saintonge et al^|2013 ). However, our model implies 
that the physical state of the cold gas varies strongly as a 
function of stellar mass along the stellar mass-SER relation. 

To understand the origin of the varying molecular frac¬ 
tions and SEE, we have to look at the galaxy properties 
that enter our model. Within our model the cold-gas and 
molecular-gas contents of a galaxy are set by a combination 
of its SER, stellar mass, and size. Most important is the 
combination of SER and scale radius, which sets the SER 
surface density distribution of the galaxy. Within our model, 
low SER surface densities result in low H 2 surface densities. 
However, to ensure the pressure of the gaseous disc is high 
enough for the cold gas to collapse and form molecules, low 
H 2 surface densities require relatively large cold-gas surface 
densities. This argument can also be reversed: an increase in 
the scale radius of the gaseous disc lowers the surface den¬ 
sity of the cold gas, naturally causing less H 2 to self-shield 
and transform into stars. We see this in the galaxies with 
low stellar masses and large discs in our sample. These are 
galaxies with the highest gas fractions, but their molecular 
fractions and star-format ion efficiencies are very low. 

An increase in galaxy scale radius also has a signifi¬ 
cant effect on the SEE of galaxies with high stellar masses 


^ We again emphasize that in this work the star-formation effi¬ 
ciency is defined as the ratio between SFR and total (HI -h H 2 ) 
cold gas mass, rather than just the molecular hydrogen mass 
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(M, > 10^° M 0 ). At fixed stellar mass galaxies with com- 
pact discs have higher SFEs (and molecular fractions) than 
galaxies with more extended discs. The lower surface densi¬ 
ties allow for a higher shielding rate of the molecular hydro¬ 
gen. At fixed stellar mass the SFE of galaxies also increases 
with increasing SER. The change in SFE cannot account for 
differences in SER of a few orders of magnitude, but it can 
lower the SER by a factor of a few. To first order it is the 
absolute amount of cold gas that largely controls the gas 
density and SER of a galaxy. The distribution of gas into an 
extended disc can act as a secondary mechanism to suppress 
the formation of molecules and lower a galaxy’s SER. 

In reality the formation of molecules does not solely de¬ 
pend on the surface density or midplane pressure of cold 
gas. Dust grains and metals are the primary catalysts and 
coolants of H 2 formation. Dust grains also shield the molec¬ 
ular hydrogen from H 2 dissociation by UV photons and 
photo-electric heating. The metallicity of the ISM is there¬ 
fore an important extra component that controls the forma¬ 
tion of molecules, especially in low-density environments. 
In addition to low densities, a low-metallicity environment 
therefore stimulates the build-up of large atomic gas reser¬ 
voirs, preventing the cold gas from condensing and forming 
stars ([Krumholz k, Dekel|2012b| [Berry et al.|2Q14| [Popping, 

Somerville fc Trager||2014| [Somerville, Popping Trager 
2015| ). We have not included this effect in our model, but 
it would strengthen the role that low surface densities and 
the partitioning of cold gas into HI and H 2 plays in the for¬ 
mation of stars. Furthermore, we have not discussed the de¬ 
struction of molecular hydrogen by radiation from an AGN. 


6 SUMMARY & CONCLUSIONS 

We applied a method to infer the total cold gas and molecu¬ 
lar gas content of a deep sample of galaxies from the CAN- 
DELS survey covering a redshift range of 0.5 < z < 3. All 
data results in this paper are available for download online 
Our main results are as follows: 

• There is an increasing trend between the inferred cold 
gas. Hi, and H 2 mass and the stellar masses of galaxies, 
which is already in place at z = 3.0 and extends down to 
stellar masses of 10® M©. The slopes of these trends are dif¬ 
ferent from the slope of the relation between SER and stellar 
mass. 

• There is no simple one-to-one mapping between inferred 
cold gas mass and molecular-hydrogen mass. Large gas reser¬ 
voirs do not necessarily lead to large H 2 reservoirs. The 
molecular fractions of cold gas increase with increasing stel¬ 
lar mass and look-back time. 

• The cold gas fraction (fgas) and relative amount of 
molecular hydrogen in galaxies (/h 2 *) decrease gradually 
with time at a relatively constant rate, independent of stellar 
mass. 

• The mean fraction of atomic hydrogen {Mhi /-f 
Mhi) in galaxies at fixed stellar mass stays remarkably con¬ 
stant with time over the entire redshift probed. 


® http://www.eso.org/~gpopping/Gergo_Poppings_Homepage/ 
Data.html 


• There is a large population of low-stellar-mass galaxies 
that are dominated by atomic gas. Only a minor fraction of 
their total gas content is molecular and can form stars. 

• The inferred SFE of galaxies increases along the relation 
between SER and stellar mass. Although part of the same 
trend, the most-massive galaxies in our sample have SFEs 
more than twice as high as lower-mass galaxies. At fixed 
stellar mass the SFE and molecular fraction of the cold gas 
increase with galaxy compactness. 

• The adopted approach can be tremendously helpful 
in understanding the impact of selection biases for much 
smaller samples of directly observed gas masses, as well 
as extending scaling relation between gas mass and other 
galaxy properties to ranges and redshift difhcult to reach 
through direct observations. 
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